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Abstract: We investigate a spinless fermion system on a one dimensional lattice interacting locally with the 
optical modes of a quantized phonon field: the Holstein model. The system is shown to have a disordered 
Q\ , ground state, for small enough coupling, at any density. This is in contrast to the non quantized phonon 

case, the static Holstein model, which at half filling has an ordered ground state for all couplings. 
>v 



§_Z Model and results. 

The hamiltonian is the sum of a free lattice fermion hamiltonian Hp, a free phonon field hamiltonian Hg 
and a local fermion phonon interaction V: 
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*I3 ■ where a? is a point on the lattice with spacing I, x = jl, j € (— Jj, ^j], j integer. Periodic boundary conditions 
are imposed at ±U (by identifying such points). The ai are creation and annihilation operators for a spinless 
fermion at x, ng = aiaZ and a+a^ + a^a+ = 5 xy . 

The field (fs represents a quantized bosonic field, corresponding to a discretized vibrating string with 
linear density a, optical frequency u and maximum wave propagation speed c. The physical meaning of <pg 
is that of deformation of the crystal cell sitting at x. The parameter m is a scale parameter fixing the bare 
fermion mass. The fermions chemical potential is fi — fiQ + v. 

The reason for writing /i in this form is that, since the early works on the theory of Fermi systems, [LW] , it 
has been realized that it is more natural to study the properties of such systems when the interaction strength 
A is varied at fixed Fermi momentum rather than at fixed fi. For the free particle system, corresponding to 
A = 0, the Fermi momentum p F is defined as the value of the momentum where the momentum distribution, 
which is the Fourier transform of the one particle reduced density matrix or equal time pair Schwinger 
function, has a discontinuity. This manifests itself in the one particle density matrix having an oscillatory 
decay Gq(x) ~ | a?| 1 sinp F |a;| andp F is obtained from \i by the relation /i = 2(1 — cos lp F ). For A ^ a Fermi 
momentum P(fi, A) can still be defined via the position of the singularity of the Fourier transform of the one 
particle reduced density matrix. In particular, for A small we shall prove that the one point reduced density 
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matrix has an oscillating leading asymptotic behaviour G{x) proportional to |x|~ 1 ~ 2?) (0(A) + sinP|cc|) for 
some T) > 0. We now define p =2(l — coslp F ) and choose v such that P(l-t, A) = p F . This defines v as a 
function of p F and A with v = when A = 0. 

It follows from our analysis (see [BGPS] and below) that P is smooth in p, A for p near any prefixed 
S (0,2) and A small enough (how small depends on the value of /xo)- Therefore we can write P = 
p F + c\(p — /xo) + biA + . . . and ci 7^ (in fact c\ = (2lsinlp F )~ 1 ). Setting P = p F then yields v = d\\ + . . . 
with d\ = b\/c\. 

We note that, according to the (formal) Luttinger theorem [L], fixing P(p,X) is equivalent to fixing the 
physical density p = p(p + v, A), i.e. P is independent of A if p is fixed: in fact P — np. The (formal) 
extension of this theorem (proved formally in [LW], [L] and formulated there as the "conservation of the 
Fermi surface volume" at constant density) to cover the present case is discussed, for completeness, in the 
appendix. 

In general the value of v is a complicated function of A, which can only be determined order by order in 
perturbation theory. There is however an exception; in fact, setting eg = (— 1) X ^ L , the unitary transformation 
% — ► £x a |\ —> —fx maps the hamiltonian with chemical potential p — 2 and v = uq into that with 
p = 2 and v = —v a and a state with density p into one with density 1 — p. Hence we see that, if v = 0, 
there must be a ground state with density |. Furthermore the hamiltonian has other symmetries; namely 
translation invariance and reflection (parity) symmetry. Thus, if we suppose that the ground state is unique 
(a property that we expect but do not prove) then by applying to it the above three symmetries we see that 
G(x) = (a^ag ) = for x/l even. Hence if one can prove, as we claim here, the existence of P such that 
G(x) is proportional to |a?| _1_2,? (0(A) + sin P|x|), it follows that P — tt/21 (otherwise, for x/l even and large 
enough, sin Pa; could be of order 1), so that the Luttinger theorem is automatically satisfied. 

The units can be fixed so that: 



h = m = I = mauj = 1 



(1.2) 



Setting b = civ 1 ,<Jq = cr, H becomes: 



2H 



J2(~ a s a z+i - 4 a 3-i + ( 2 - m)4% ) 



(1.3) 



with the half filled band case corresponding to p, = 2, v = 0. 



§i.i The Static Holstein Model. 

If (To = +00, 6 = the model (1.3) is the static Holstein model with "magnetic field" equal to A/4, [LM]. 
The ground state problem is now equivalent to the computation of the fermionic energy E(X, v. <p) in the 
presence of the external field ipg and then finding the field ip minimizing: 

£7(A,i/ ;¥ >) + ±5>| (1.4) 

X 

Calling £(A, v) the minimum value, the corresponding density is obtained by setting p — \ = —2d v £{\,v) 
which, if p is given a priori , is an equation determining v as a function of A, p. 

The case p = 2,v = 0, p = ^, the so called "half filled band" case, has been solved by [LM], following 
the methods used by [KL] for the case in which the field ipg can only take values ±1 (the Falikov-Kimball 
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model, see also [MM]). It is shown in [LM] that the minimizing field is ip$ = (— l) x f where / is a suitable 
constant, which can be easily computed by remarking that E(X, v; (—l) x f) can be evaluated by the Bloch 
waves techniques. 

One finds that A ^ implies / ^= 0. This means that, if A 0, there are two non translation invariant 
ground states, in which the field (ps has a periodic structure with period 2. This is interpreted by saying 
that at half filling the atoms of the system acquire a crystalline ordering with a period 2, at any non zero 
coupling strength; i.e. the Peierls instability. They, thus, behave as if they were non interacting particles 
immersed in a periodic potential with period 2. 

It follows from this that the one particle energies are split into two bands, the first corresponding to the 
momenta \k\ < |, and the second to the larger k's. The two bands are separated by a gap k(X) and, as a 
consequence, the fermions are in a filled band state (because the Fermi momentum is p F = ^ for p = 1/2, 
by the above symmetry). Consequently the two point equal time Schwinger function for the fermions, 
G{x) = lim t „> - G(x, t), decays at large separation as, 

\G(x)\ ~ e - K( ^ (1.5) 

where k(X) is the energy gap at momentum p F = ir/2. Note that the energy gap k(X) a ^ > 0i and that at 
A = one has instead: G(x) = Gq(x) = —(nx)^ 1 smp F x. 



§1.2 The Dynamic Holstein Model. 

The above picture holds for all values of A ^ at half filling (p — \, v = 0) when the phonon field is 
treated classically, by setting a^ 2 = in (1.3) from the beginning. We shall show that the results in [BGPS] 
imply that, if Co < oo, i.e. if quantum effects arc not neglected, then for any density p e (0, 1), and A small 
enough (depending on tro, b), the decay of G(x, 0) is given by: 

G ( x >°)^^m (L6) 

with 77(A) analytic in A 2 ; 77(A) = aA 4 + (3(A 6 ) if b ^ 0, or 77(A) = a' A 6 + (9(A 8 ) if b = 0, with a, a' ^ 0. 

This is clearly incompatible with a periodic minimizing state of the field ip x , showing that the Holstein 
model ground state for a spinless fermion system is "disordered" at small coupling, i.e. there is no long range 
order in < ipgipg > or < ngn^j > as \x — y\ — > 00. The situation could change at large coupling, where the 
ground state could again be ordered: but this is outside the domain of applicability of our techniques. The 
maximum value of A, |A| < A m ((T ), for which we can prove (1.6) goes to zero as a — > 00. 

The problem of fermions with spin is, of course, much more interesting. By the same argument given 
below this case can also be reduced to the problem of fermions with spin, interacting with a short range 
potential. However the fermionic interaction is attractive: and while this has no consequences in the spinless 
case it has profound consequences in the case with spin. Unfortunately the spinning case is not understood 
in the sense needed here to draw any conclusion, sec [BM]. 



3 



$2 Reduction to a continuum problem. 



When A = v = the phonon and fermion fields are independent. Their respective Schwinger functions 
(imaginary time Green functions, see [FW] for definitions) can be computed via the Wick rule from the two 
point functions Sa(x) and Qa{x), where A=[0, /?] x (—7, 7] and x = (t,x) 6 A: 



1 a —ikot—ikx 



ikot— ikx 



(2.1) 



^(*)=*7 E E 
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P L e^=- le ^ =irkl<7r -^o + e F (k) 



where k = (fen, k), and: 



es(fe) = 2(1 — cos k), ep(k) = (cosp F — cos k) (2.2) 

The summation rule over fco is to take the limit o£ | fe 1 < JV as ^ ~~ * °°> an0 - we su PP ose that p F = ^(n^ + i), 

with rif integer, so that ep-(fc) 7^ for all k. 

An application of Trotter's formula allows us, as usual (sec [BG1]), to write an expression for the ground 
state infinite volume Schwinger functions of the interacting fermions, G(x), in terms of the gaussian integral 
Ps(d$) with propagator Sa in (2.1) and of the grassmannian integral Pp(rfV') with propagator g A in (2.1); 
for example: 

J P B (d^)P F (d^)e- v ^^+ 
j PB{d<S>)P F {d4>)e- v ^^) 

V A (*,il>)=[ dt Yl (-^ + A$x)^C (2-3) 

Se(-L/2,L/2] 

where the fields $ satisfy periodic boundary conditions in x and t, while the tp fields (which are grassmannian 
variables) satisfy periodic boundary conditions in x and are antiperiodic in t. 1 The term (A 2 /2)S\(0)ip+'ip~ 
in the definition of Va is produced by the shift of the field & x which makes it possible to write the measure 
Ps(rf$) as a measure with zero mean. 

The propagator S(x) obtained from S A (x) in the limits (3 — > exo, L — > 00 is exponentially decreasing; in 
fact, if \t\ < (3/2 and \x\ < L/2 , one can easily see that 2 : 

\Sa(x)\ < < ^- e - Klltla ° 1 e- K2 ( b)lsl (2.4) 
where Ki < 1 is independent of b, and: 



fOQog^ 1 ) for&^O (0(1) for 6^0 

^ W- \0(6- 1 ) for&^oo ' W \ ©(r 1 log ft) for&^oo 

In particular, for 6 = it is, in the limit (3 — > 00, S(a;) = (2cr ) _1 e _ l t l' T f) 



(2.5) 



1 One should not confuse the fields <I> (which are real valued fields) with the operators ip x ; nor should one confuse the 
grassmannian variables ip with the operators a x . 
Recall that the boundary conditions arc periodic. 
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The integral over the bosonic field $ is a gaussian integral, which can be performed explicitly, yielding: 



and, denoting J A dx- = Exg(-l/2,l/2] Jo 



V A (il>) = -v J^dxi>+il) x - y J^dxdyS(x-y)ip+ip x tp+ ip y (2.7) 

We see that the problem has the same formal structure as that of a fermion system on the continuous 
line interacting via a short range potential, considered in [BGPS], see [Gl] for a summary. The attractivity 
is due to the positive definitness of S(x), see (2.1). 

Let us point out the main differences: 

(1) In formula (2.3) the propagator g(x) for the grassmannian integral, defined in (2.1) and (2.2), has a 
different "dispersion relation"; in[BGPS] the dispersion relation is ep(k) = \{k 2 — p 2 F ), while in the present 
case eF(fc) = {cosp F — cosfc). The difference is due to the fact that in the Holstcin model the fermions are 
on a lattice. 

(2) The fermion potential (2.7) is non local in time, while in space it may even have zero range (in the case 
6 = 0, sec (2.5)). In the continuum problem considered in [BGPS] the potential had the form (2.7) with 
v(x — y)S(t — t') replacing S(x — y), with x = (t, x), y = (f, y). 

These changes are of no consequence: in fact what was really used in [BGPS] was that ^-(Pf) > which 
is true in the present case as well, provided < p F < n, as we suppose ("positive density smaller than close 
packing"). 

Therefore we can perform the decomposition of the propagator in the same way as in [BGPS], see (13) in 
[Gl], by writing: 

1 1 - e - fe o+ e F(fe) 2 e -k 2 -e F (k) 2 



-iko + ep{k) —ik + ep(k) —ik ~eF(k) 



(2.E 



(2.8) and (2.1) generate a decomposition of the propagator into a sum of an ultraviolet part, g u , v \x)= 
g(>°)(x), and of an infrared part, gi. r .{x)= g(-^{x). The decomposition (2.8) allows us to represent ipf as 
sums of two independent grassmannian variables, that we denote v4 >0 ^ and ipi^ ' 1 with respective propagators 
0<>°) and g^°K 

The integration over ^( >0 ) can be controlled by perturbation expansions, as in [BGPS], and we shall 
integrate over it. The remaining integral for the evaluation of the partition function, i.e. for the denominator 
of (2.6), then becomes: 

" P F (d^)e- v(a) ^- 0) ) (2.9) 



with: 



/ 

y(°)(V>) =- y J dx dyS(x-y)iP+iP+^-ip- + 

-A 2 J dxdyS(x-y)g ( > >a ')(y-x)^i>y + 

-(P + AV >0) (0)) I dxi,+ij-+jr fdx 1 ...dx 

J n=l 
■ W 2 „(X 1 . . . X 2n )^t 1 ■ ■ ■ tt„ +1 • ■ • C 2n 



(2.10) 
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where the kernels W2 n {xi ■ ■ ■) are analytic in (A 2 ,^) = r for \r\ < e (for some e > 0) and verify the short 
range property: 

J dx 1 ...dx 2n \W 2n (x 1 ...x 2n )\e Kd ^-> x ^ < \A\(D\r\r^ 2 < n -V (2.11) 

where k = \ min(Ki<jQ K2(b)), see (2.4), and d{x\, . . .) is the length of the shortest path connecting all the 
points xi, . . . (regarded as points on the torus A). 

The main result achieved by (2.10) is that we have "disposed" of the ultaviolet part of the problem (in the 
evaluation of the partition function) and we can say that the problem is reduced to an essentially identical 
one with a purely infrared propagator and a new which, to lowest order in the couplings r = (A 2 , v) 
has the same form as the original one, as far as the quartic part is concerned, and a slightly different (non 
local) quadratic part; plus "higher order terms" of every degree in the ip fields. All the terms in (2.10) are 
well defined, and have a convergent power series in r = (A 2 , v) if r is small enough. 

The proof of the above statements is simpler than the corresponding one in [BGPS], see theorem 1, because 
in the present case there is a natural ultraviolet cut off, at least in the x direction, so that no multi-scale 
decomposition of g u .v.{x) is needed. 

The non locality of the quartic part in V^ ^ is not important even for the infrared problem as it does not 
affect the notion of relevant or marginal and irrelevant operators, which is the notion on which the infrared 
integration is based when performed via the renormalization group method in [BGPS]. Hence we reach the 
same conclusions about the partition function and about the Schwinger functions (whose analysis can be 
performed once the partition function can be estimated in detail, see [BGPS], §5 and §6). The beta function 
is essentially (as shown in §7 of [BGPS]) and this allows us to draw the "same" conclusions as in [BGPS]. 

In particular one can prove the anomalous asymptotic behaviour (1.6) of the two point Schwinger function 
with a coefficient 77(A) which is analytic in A 2 and in general is of order A 4 . However, if b = 0, one can see 
by an explicit calculation that the leading term in the expansion of 77(A) vanishes, while the following one 
seems different from zero. 

Another interesting observation is that the formal "Luttinger Theorem" [LW] is valid in our case: it states 
that the density p is a function only of p F , that is p — p F /n for any A: see appendix for a discussion of this 
formal result (for which it would be nice to have a rigorus proof). 

In the next section we give some details on the recursion procedure that we follow to obtain the estimates, 
and on how we define the relevant operators: this will follow closely [BGPS] but should be useful to the 
readers who wish to see where one is going, before plunging themselves into the analysis of our estimates. 
But neither the estimates nor the proof of the vanishing of the beta function will be reproduced here: they 
are identical to the corresponding proofs in [BGPS]. 



§5 The recursive evaluation of the partition function. 

The following is the "standard" anomalous dimension renormalization procedure. It is illustrated in quite 
a general context, in which one is given a priori an arbitrary sequence of "wave function renormalization" 
constants: Zq=1, Z_i, Z_2, 

The infrared propagator g(-°\x) is decomposed as: 



9 ( ' J (2^) 2 - iko + e F (k) J (2^) 2 - iko + e F {k) 1 ' 

where t h (k) = cxp(-2- 2,l (fcg + e F {k) 2 )) and T (k) = t {k) - 

Eq. (3.1) allows us to represent ^>(-°) as tp^ + ipt-^ 1 * 1 where t/^- -1 ) has propagator Z ( ^" 1 i_i(/c)(— iko + 
epik))- 1 andi/»(°) has a propagator given by Z^T^k^-ikQ + e F (fc))- 1 . Let P Zo (dtp^) and P Zo (d^- 1 ^) 
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denote the corresponding integrations, the grassmannian integral in (2.9) thus becomes: 

? Zo (#(°))Pz (#^- 1) ) (3.2) 
Using the sequence Zj one can then write the identity: 



{Z Q tz{ + Z-i - Zo)- 1 = ZZ{t- 2 + ZZ\T { 1} (3.3) 



defining T ( 1} (fc). 



Hence if ip( ^ is a grassmannian field with propagator: 

1 



Z-i -ik Q + e F (k) 



(3.4) 



and ip^- 2 ' is a grassmannian field with propagator Z-it-2(k)(—iko + ep{k)) 1 , independent of ip( ^ and 
of ip(°\ then the grassmannian integral (3.2) becomes (formally), up to normalization constants: 



e 



(Z_i-Z ) / 4>i-~ 1>+ (dt+eiid?))^- 1 



(3.5) 



as is easily checked with some algebra, if ip(- ^ = t)A ^ + 2 '. 
By iteration: 

w ( - 0) )= n? Zl w (j) ) ■p Zh (d^- h - 1) y 



h+1 ,. 

cxp^(^_ 1 - Zj) / ^- l)+ (d t + e{id s ))^~V- dx 

3=0 J 



(3.6) 



and?A^ ) = ^ (0) +^ ( " 1) + ... + ^ ( ' l) +^ ( - /l " 1) , withV> (i) having propagator ZJ 1 ^ {k) {-ik + e F (fc))- 1 ; 
Z^gU) w ith: 



r (j) (fc) = ij - + (1 - TTZ7T 

if is defined by: Zj=(l + Zj)Zj + i. 

The integration (2.9) can therefore be performed recursively by setting: 



(3.7) 



e- vW (^ ( -")= f f[ Pz h ,(d^ 

J h' = h+l 



>=h+i (3i 



so that: 



| e -v(°) W <^) pW (<o )) = | e -^)(^<^)p ZoW (<o) ) 
= Jp Zh (d*P {h) )Pz h (d^ <h) )e- 



-v {h) (^zZip { - h) ) 



(3.9) 



where = + V^ -1 ^ 

The (independent) "fluctuation fields" tp^ have propagators with good scaling properties if \zh\ < \- 
they can be represented via quasi particle fields (see below) with propagators bounded, unformly in h, by 
2 h ^(2 h x) for some function 7(2) which decays exponentially fast as x — > 00. 
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The idea is to select the sequence Zh so that \zh\ < \ and so that the potential does not contain 
certain terms which would otherwise be difficult to control. 

To understand the choice of Zh one has to define the relevant and marginal terms in V^ h \ Such notions 
are not naturally defined from the considered as functions of the particle fields ip(- h \ They are very 
natural if V (h) is regarded as a function of certain auxiliary fields that we call the quasi particle fields. 

The infrared propagators g(- h \x) can be decomposed as: 

9 { - h \x)= E e- ipF " s g^ h \x,uj) 

u5=±l 

(3.10) 

<<h) (x ~ =e i PF ** r dh r dk t h (k) X (2- h zk) c _ ikx 

J-oo 2tt J_ n 2n - 2 fc + eF (fc) 

where x(r) may be chosen as 7r -1 / 2 e~ s2 ds, so that x is a smooth version of the step function, with the 
property that x(r) + x(~ r ) = l- 

It is easy to check that, in the infinite (space-time) volume limit, g^- h \x,Q) is essentially scale indepen- 
dent, i.e. h independent, as h — > — co. In fact: 

g^ h \x,uj) ^fc^_oo 2 h g(2 h x,uj) 

dkt{kl + (,inp F fkf)^ lkx (3- 11 ) 



^ _ r°° dk r°° dk 
J-oo 27r J_ 00 2tt 



-iko + Ldksmp F 



-e 



The decomposition (3.10) of the covariance allows to perform a related decomposition of the fields ip^~ h>> : 



i>t ( ~ h) = E e^HT h) (3-12) 

<3=±1 



The fields ^~~ h \ which will be called the quasi particle fields, are independent fields with propagator 
g(< h )(x, w); moreover they are antiperiodic also in the x variable, so that we can write: 

v4= E e±ikx ^ ( 3 - 13 ) 

e ik /3 = _ le ikL = _ 1 

One can think that the ip^-^ have a distribution which, "up to scaling", is h independent. This means 

that the distribution of ip(- h ^ is the same as that of 2 h ^ 2 tp 2 h x.i2 where ?p x ^ has a propagator g(x, uj)5qs', see 
(3.11) (up to corrections vanishing fast as h — * —oo). The fluctuation fields can also be expressed in terms of 
quasi particle fields in the same way: their propagators have the good scaling properties mentioned above. 

Thus if we think of the infrared problem as that of integrating over the quasi particle fields tp^ h \ we 
have a natural way of introducing the notions of relevant, marginal and irrelevant operators. Precisely we 
define the relevant operator to be: 

p 2 = E / ^^ i( *-*' )PFS dx = J ^~dx (3.14) 

and the marginal operators as: 

F 4 =J TP+ + ^-_ip- + dx (3.15) 
F 2 , T = E / e 4( ^ >Fi Vi(d t + iu' smp F V s ,)i,- a , dx = J + e F {id s ))i>- dx 

F2.a = E / S*- '** 3 ^ iZ'sinp F V^- c3 ,dx = I ^+e F (id s )^- dx 
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where is a suitable operator, which acts by multiplication by M$(k) on the Fourier transforms, with 

M a (k) = -2ismk/2 Sin{PF+ " k/2) ~ g n -zfc (3.16) 

Note that, while the second degree operators can be expressed easily in terms of the particle fields, the same 
is not true for the F4. 

The "usual" power counting attributes a size to the above operators evaluated by extending the integral 
over a box of size 2~ h , i.e. of volume 2~ 2h , and by attributing to each field a size 

2 h/2 

as suggested by the 

scaling properties discussed above; furthermore each derivative contributes to the size a factor 2 h . 

Hence the conventional power counting attributes to F2 a size that is evaluated as 2~ 2h (2 h / 2 ) 2 = 2~ h 
(hence F 2 is relevant). The size of F 4 is 2- 2h (2 h / 2 ) 4 = 1 and F 2 , r ,F 2i<J have size 2- 2h 2 h {2 h ' 2 ) 2 = 1 (hence 
they are marginal). All the other local operators are irrelevant (i.e. they have sizes 2 h ' 2 or less). 

Note that the size of the various operators is clear if they are regarded as functions of the quasi particle 
fields. 

Given V^(ip), we must identify the relevant and marginal parts of V^ h K This is done by introducing the 
localization operator C: it is a linear projection operator which is unless acting on a fourth degree or on 
a second degree monomial in the fields. Imagining that is expressed in terms of quasi particle fields, 
hence as a sum (or integral) of monomials in the quasi particle fields, the action of C on the fourth degree 
monomials is described by: 



1 2 

Cib+ . xh + . th- , V" - = - yV> + - ip~ - (3.17) 



while the action of C on the second degree monomials is: 

Wta^ySi = ^tsMxS 2 + Mw> - XQ)d Xo i)- a ^- a h L {y-x)V a ^- w ^ (3.18) 

where 

h M (s -t) = il[ e Ws-*)/ M _ e -i*(«-t)/M] (3 19) 

is an antiperiodic approximation of (s — t), which converges to it as M — > oo. 

The operator T>£ differs from the analogous operator defined in [BGPS], where x was a continuum variable. 
However, it is easy to see that one has to change the bounds obtained in [BGPS] only in minor points, without 
affecting the final results. In (3.18) we have also taken explicitely into account the boundary conditions, 
while in [BGPS] this problem was neglected. However, as explicitely shown in the analysis of the spinning 
case in [BM], this approximation does not play any relevant role. 



Appendix 

In this appendix we give a heuristic derivation of the so called Luttingcr Theorem, following the analysis 
presented in [LW]. We believe that the techniques used in [BGPS] to study the Schwinger functions could 
be used also to prove rigorously the Luttinger Theorem, but this has not been done yet. 

We consider a system in finite volume L at zero temperature. If we define the hamiltonian H so that the 
ground state |0) has eigenvalue 0, the two point Schwinger function is defined in the following way: 

G(x,t) - X (t > 0)(0\aZe- tH a+\0) - X (t < 0)(0\a+e+ tH aZ\0) (Al) 
where x(A) is the characteristic function of the set of points where the condition A holds. 
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Hence its Fourier tranform G(k, ho) is given by: 



G{k,k )= ( dt\e lkot {0\aZe~ tH ai\0) -e- lkot {Q\a-±e- tH aZ\Q) 

= Y / dt le-^-^ma^n)^ - e-*^^\(0\a+\n)\ 2 

£ — ' /. [ k k 

n:E n >0 J() 



(A2) 



y r |(0|ag|n)| 2 | |(0|4|n)| 2 - 
I -ik + E n -iko - E n 

where E n , n = 0, . . ., are the eigenvalues of H and \n) the corresponding eigenstates. The ground state gives 
no contribution to the sum in (A2) (so allowing the integration over t) because of the conservation of the 
total momentum, which implies that (0|ai |0) = 0, when k ^ 0, and because of a cancellation between the 

two different terms in the r.h.s. of (A2), when k = 0. 

The function G(fc), k — {k, fco), is often written also in the following form: 

G(k) = L (A3) 

-ik + e F (k) + E(fc) 

which defines the self-energy function S(fc). As is well known, [LW], S(fc) can be expressed in perturbation 
theory as the sum of all connected graphs with two external lines, which are irreducible, that is which do 
not become disconnected by erasing any internal line. 

By using (Al), one sees immediately that the density p is given by: 

p = - lim G(0, t) = ^ Y lim / ^ (A4) 



fe 



Let us write 

^-^ = £ log[C - e F (k) - m -iC)] + ^ (g '- <)/9C (A5) 

by choosing the branch of log z so that logz = log \ z\ + i&rgz, < argz < 27T. 

It was argued in [LW], if we insert the r.h.s of (A5) in the r.h.s. of (A4), the second term gives no 
contribution, that is: 

iVlim / ^e^ mwdko i w ^ s(fc) g|w = Q {M) 

L^e^o+J 2ir ik Q -e F {k)~X(k) L^J 2n dk 

We recall briefly the arguments given in [LW] to justify (A6). Indeed, our problem is not explicitly 
considered in [LW] , where the interaction between the fermions is local in time; hence only the case b = is 
covered by [LW]. However, their arguments extend in a trivial way to our general case. 

We start from the observation (not easy to prove rigorously) that can be thought as the sum of all 
graphs (connected and irreducible) with two external lines, such that the internal lines carry the complete 
propagator G(k) and the following condition is satisfied: there is no proper subgraph with two external 
lines. We shall call such graphs, as usual, skeleton graphs and say that a graph is of order n, if it contains n 
four- fermions interacting terms (hence it has 2n vertices). 

Let Ti n (k) be the sum of all skeleton graphs of order n and let Y n be the sum of all vacuum (that is with 
no external lines) skeleton graphs of order n, which can be obtained by closing the two external lines of a 
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graph contributing to S(fc) with a propagator G(k). Since each graph contributing to Y n can be obtained 
in 2n different ways from a graph contributing to we have: 

k 

Moreover, given a graph Q contributing to Y n , its value can be written in the following way: 

Y° = J dkW . . . dfc (2n) d<Z (1) . . . dq (n) G(k ( V) ■ ■ ■ G(k (2n ^)- 

■ S(q (1) ) . . . S(q (n) ) Y[ 5(fc (si) - fc (n) + q^) 
i=l 

where J dk=L~ x J rffco/(27r) and fcW, . . . , fc( 2n ) are the (space-time) momenta of the 2n fermion lines, 
. . . are the momenta of the n phonon lines and, finally, k^ Si ^ and k^ are the momenta of the 
fermion entering and leaving the vertex i, respectively, while q^ is the momentum of the phonon propagator 
entering vertex i. Moreover, if we consider the sum of all the quantities that we obtain by substituting, in 
the r.h.s. of (A8), one fermion propagator by its derivative with respect to fc , we get: 

2n 



V / dk^ . . . dk^dqW . . . dq^G(k^) ■ ■ ■ -^y(fcW) • • • G(k^)- 
-=i J <9/cq 

2n 

S(q W ) . . . 5(<z (n) ) Y[ 8{k M - k^ + q (bi) ) = 
i=l 

2n „ 

= ± E/ dk(1) ■■■dk {2n) dq^ ...dq^G{k^)---G{k^)- 



(A9) 



d d 2 " 



• s ^ • • • + 3) ) n ^ - fc(n) + ^ (6i) ) = ° 

dk K 3 ' dky> j=1 

If we now sum the l.h.s. of (A9) over (/, we get, by an argument similar to that used in order to prove (A7), 
that: 

/BG 
dkZ(k) — (k) = (A10) 
oko 

so that (A6) is proved. 

The formal identity (A6) implies that p = j; SfePfc' 

p % = £ hm y ^ ^e^- log[C - ej ,(fc) - S(fc, -iC)] (AH) 

Moreover, (A2) implies that the integrand in the r.h.s. of (All) is analytic in all the complex £ plane, except 
on the real axis, where there are branch points at ( — ±E n , n > 0. Therefore, we can deform the integration 
contour into the contour C of Fig. 1 
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Fig. 1: Integration Contour for (All). 

and then we can integrate by parts. We get: 

= j. log[-e f (fc) - S(fc, -V)} - log[-e F (fc) - g)] 
P " ^ ^ ' (412) 



- lim £ / i^e* log[C - e F (k) - E(fc, -iC)] 



It is easy to see that the second term in (A12) vanishes, if the argument of the logarithm has a negative real 
part for Re£ — > — oo, by recalling that we have chosen the branch of log z with the cut along the positive 
real axis. 

We want to show that this is certainly true, if |(0|oi|n)| — > sufficiently fast, as E n — > oo. We note that 
Re [£ — ep(k) — £(fc, —iQ] has the same sign as — ReG(fc, — iQ; moreover, by (A2), if £ = u + iS: 

/+oo u — E 

^ (J k(dE) [u _ E)2 + 52 (413) 

where <j^{dE) is a probability measure, as it is easy to check. Moreover, if u < —uq, with uq large enough, 
there is an interval [E 1 ,E 2 ], such that \u — E\> \u\/2 for E e [Ei,E 2 ] and Jj B1 E ^ a^(dE) = 1/2; hence it 
is easy to show that: 

-ReG(k,-i() < (26)- 1 fjkidE) -J^ (414) 
which immediately implies that — ReG(fc, is definitely negative for u — > — oo, if 

lim u [ uAdE) = (415) 

Since, by (A2), the function E(fc, —iQ is real at £ = 0, it follows from (A12) that: 

P= | X>M*) + £(0, fc) <0] (416) 

fc 

Equation (A16) implies that, in the limit of infinite volume, there is a relation, independent of the strength 
of the interaction A, between the density and the Fermi momentum p F , defined as the value of k, such that 
ep(k) + E(fc, 0) = (that is the value of k where the interacting propagator is singular for k = 0). For 
A = we have p — p F /n, hence this relation has to be valid for any A. 
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